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1 Introduction 



The problem of large -order behavior of the perturbative theory attracted 
a renewed attention recently. One of the motivating factors is a common 
wisdom that the corresponding asymptotic behavior is related somehow to a 
very deep physics. This is the area where perturbative and nonperturbative 
physics strongly interfere. An understanding of this interplay may shed the 
light on the nature of the nonperturbative vacuum structure in general and 
the origin of vacuum condensates in particular. 

With these general remarks in mind we would like to analyze these prob- 
lems in solvable two-dimensional QCD{N oo), We would like to 
test assumptions, hypothesis and interpretations, made in 4-dimensional field 
theory, within toy two-dimensional QCD2{N — >■ oo), where we expect con- 
finement and many other properties inherent to real QCD. Additionally, we 
extend the analysis to QCD with adjoint matter [^0. As is known, in this 
theory, the pair creation is not suppressed even in the large Nc limit, and 
thus, this model can mimic an exponentially growing density of states with 
large mass p{m) ~ expm. In this case no exact solution is available, but we 
argue that general methods such as dispersion relations, duality and unitar- 
ity can provide all information we need about spectrum for the calculation 
of large order behavior. 

Why are we so conscious about the large order behavior? We see at least 
a few theoretical and phenomenological reasons for that. Let us start from 
the pure theoretical reasons. One may think that the crucial question in this 
case is whether the perturbative series is Borel summable or not. 

Contrary to the common belief, we do not think that the issue of Borel 
summability (or its loss ) is the fundamental point. In particular, let us 
mention an example of the principal chiral field theory at large [|r^. In 
this case, the explicit solution as well as the coefficients of the perturbative 
expansion can be calculated. These coefficients grow factorially with the 
order and the series is non-Borel summable, but nevertheless, the physical 
observables are perfectly exist for any finite coupling constant. The exact 
result can be recovered by special prescription which uses a non-trivial pro- 
cedure of analytical continuation (which might be a good example for other 
asymptotically free theories). 

The second important theoretical issue can be formulated as follows: Be- 
cause of dimensionality of the coupling constant in QCD2 the perturbative 
expansion ~ Z](5'^)"c„ and the operator product expansion (OPE) for a cor- 
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relation function ~ S(q2)"'c„ are one and the same expansion. From this 
simple observation we learn that the OPE is an asymptotic series. Thus, 
many interesting questions arise: 

a) What kind of vacuum condensates are responsible for such a behavior? 

b) Do we extract the actual condensates from the OPE or only effective 
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ones? 

c) What kind of vacuum configurations are responsible for such n! growth? 

d) Do these configurations saturate the vacuum condensates? 

What is more important, sometimes these questions (and many others) 
can be answered. We expect, will argue, that the analogous phenomena 
might occur in real four dimensional QCD, thus these questions have not 
only pure academic interest. 

Phenomenologically, there are issues which are even more interesting and 
give much more freedom for speculation. First of all, let us recall that the 
reason for interest in the large order behavior is related to the factorial n\ ~ 

growth of the perturbative coefficients. This growth can be considered 
via the dispersion relation and is commonly interpreted as a reflection of the 
divergence of the multiparticle cross section with large number of particles 
and energy ~ n (see discussions in The naive interpretation would be 

the violation of unitarity 0. We will show however that in QCD2 while we 
have a factorial growth of coefficients, this growth has nothing to do with 
multiparticle production since at large N the pair creation is suppressed by 
factor Rather, this growth is related to the highly excited two-particle 
meson states. Another phenomenological issue looks mysterious: vacuum 
condensates extracted from the spectrum might be quite different from the 
actual magnitude of condensates. 

2 t'Hooft model 

Let us start from the analysis of two dimensional QCD with fermions in 
the fundamental representation -the 't Hooft model. It is completely solvable 
in the limit where the number of colors — > 00 The Bethe-Salpeter 
equation for mesonic bound states was solved in |I| yielding a spectrum 
whose states lie asymptotically on a single "Regge trajectory". We want to 
point out that many general questions in this model can be answered without 
solving an equation, but using such powerful methods as dispersion relations, 
duality and unitarity. In particular, in the weak coupling regime, 

fif^A^ ~ const. N ^ 00, nig^ q —= (1) 

^ \/iV ^ ' 

the chiral {t^)'^) and gluon condensates (G^j^) can be calculated exactly., see 
below. Additionally, few low-energy theorems can be tested and obtained 
result imply that there are no other states in addition to those found by 't 
Hooft. In other words, the dispersion and duality relations would indicate 
missing states. 

^In the physical theory, the unitarity is preserved, of course. The physical question is: 
what can stop this growth? 
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Here, the entire spectrum is discrete and is classified by the integer n. 
The model we shall consider consists of quark in fundamental representation 
interacting via an SU {N) color gauge group. We follow the notation of ref.[^ 
and present the 't Hooft equation |T[] in the following form 

where symbol P notes as the principal value of the integral, and < x < 1 
is the fraction of the total momentum of the bound state carried by quark q 
with mass m^. The quantity ml = is the basic mass scale in the theory 
and the index n classifies the ordering number of the bound states \n,p) with 
total momentum p^. 

The same wave function can be expressed in terms of the following matrix 
element 0: 

M^) = /f/ dy+e-^'^^'-'''^'-{0\q{-y)q{y)\n,p)\y^=o. (3) 

Let us note that the matrix element on the right is written in the light cone 
gauge y4_ = 0; to restore the manifest gauge invariance one can insert the 
standard exponential factor S ^-''■y+ into the formula (|^). 

Let us review some important properties of equation @. The entire spec- 
trum is discrete and classified by the integer number n. The wave functions 
4>n{x) are orthogonal, complete and obey the following boundary conditions 
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Tfi 

(pJx) ^ \x(l - X)f , X^O, X^l, 7T (3 COt( 7f l3) = 1 1. (4) 

For large n the spectrum is linear 

~ n'^mln, 0„(x) ^ \/2sin(7mx) (5) 

and does not depend on mass of the quark. More importantly, in the chiral 
limit {rriq — > 0) the lowest level (we call it vr meson) tends to zero (m^ ~ 
rriq) and one could expect a nonzero magnitude for the chiral condensate. 
Thus, we have come to the very important connection between spectrum 
and vacuum structure. 

As the vacuum of the model is a very important issue for the following 
analysis, we would like to recall some results with an explanation of the 
general methods which have been used to derive them. 

We define the chiral condensate in the current algebra terms as follows: 

= lim ^ / d^xe'P^d^{0\T{q^,j,q{x), g75g(0)}|0) = (6) 
2t{0\qq\0) + 2m, • {0\T{qt^,q{x), toWIIO). 
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As we have already mentioned, the only states of 't Hooft's solution are 
the quark-antiquark bound states, thus they must saturate the dispersion 
relation. Upon inserting this complete set of mesons to the (0) one obtains: 



{0\qq\0) 



-m, 



E 



(7) 



where /„ is defined in terms of the following matrix elements 



{0\qi-f5q\n) 



-7rmo/„, /nvrmo 



x(l — x) 



dx 



(8) 



In the chiral limit the only state which can contribute to the formula (|^) is 
the TT meson. Its matrix element can be calculated exactly and we end up 
with the following expression for the chiral condensate in the — >■ limit 
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(9) 



The result is confirmed by numerical [r^|,[T^ and independent analytical 



calculations |jT5|. Moreover, the method has been generalized for the nonzero 
quark mass and the corresponding explicit formula for the chiral condensate 
(qq) with arbitrary rriq has been obtained . 



As was expected, we find that (0|gg|0) 

if we put nig 



N. Besides that, as we have 
from the very beginning, then 
This corresponds to the different regime when rriq g ^ 



already noticed in [|T2| 
(0|gg|0) = 0. 

when nonplanar diagrams come into the game as we will discuss 
later. The last remark is the observation that the entire nonzero answer for 
the condensate comes from the infrared region of the integration in eq.(^: 
X ~ 0, X ~ 1 which corresponds to the situation when one of the quarks 
carries all the momentum and the second one is at rest. 

The sum (|^) can be calculated exactly for arbitrary rrig The crucial 
point is that for arbitrary nig the nonzero contribution comes from the highly 
excited states {n ^ 1) only. The properties of these states are well-known: 



J n 



TT TTIq ■ n 



n > 1, 



(10) 



and thus the sum (|^) can be explicitly evaluated with the result [|T^ 



m N 1 
(0|gg|0) = ^{log(7ra) - 1 - 7^ + (1 - -)[/(«) - - log 4]}, (11) 



where a = — •je = 0.5772.. is Euler's constant and 

dy 1 



Ha] 



sinh y cosh y 



y"^ cothy — 1) + 1] 
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This result is exact for large A'^ and arbitrary quark mass within the 't Hooft 
regime, i.e. ^ (7 ~ (|T]). In the limit a — » 00, it reduces to the 

eq.(^)as it should. 

The last condition {rriq ^ g ) which has to be satisfied for the 't Hooft 
solution to be valid, requires some additional explanation. Roughly speaking, 
nonplanar diagrams may contain a factor ~ which at rUq = blows up 
and the theory changes completely. The concept of the proof that there exists 
a factor ~ in nonplanar diagrams is the following. 

Let us consider the correlation function for p — > 

i J ciV(0|T{gg(a;), gg(0)}|0) = P(p') (12) 

The 't Hooft solution suggests that only planar graphs are taken into account 
and, consequently, the spectral density contains only contributions from one 
meson states for which Ppianar ~ N . At the same time in the chiral limit, 
we can calculate the two-pion contribution exactly! This contribution is not 
accounted for in deriving (0). Of course, the two-pion contribution is 
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suppressed by a factor 1/A^. However, it contains a term ~ ^ which tends 
to infinity for — 0. The presence of the factor ~ in nonplanar 

diagrams leads to the aforementioned constraint on m^. 

Now let us explicitly demonstrate the existence of the term ~ for the 
two-pion contribution. In order to do so, we write down a dispersion relation 
for P: 

P(0) = - / —ImP{s), (13) 

71 J Ami ■S 

where ImP{s) is the physical spectral density. The vrvr contribution is fixed 
uniquely by (|^) because of the special role of pions [|T^: 



V3 ^ 6 Js(s-Aml) 

P^^(O) = 



^2^2 „oo mpTT^ 1 



6 ^4m| sJsis-Aml) 12m2 ruq 



It is clear that the only cause for a singular ~ l/rrig behavior is the finiteness 
of the pion matrix elements at zero momentum. At the same time this 
contribution does not contain the large factor which accompanies a one 
meson contribution to the same correlator. To suppress these nonplanar 
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diagrams we require N ^ Thus, we expect that some kind of phase 
transition may occur in the region rriq ~ g, which would couse a complete 
restructuring of the theory. 

The last subject we would like to discuss in this section is the strict 



Coleman theorem [17| which states that a continuous symmetry cannot be 
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broken spontaneously in a two dimensional theory. As we discussed earlier 
we expect that as in the SU{N — > oo) Thirring model (where the chiral 



symmetry is "almost" spontaneously broken |jT8[), the Berezinski-Kosterlitz- 
Thouless (BKT ) effect operates in regime (|ip. This fact also confirms 
the 't Hooft spectrum: states with opposite P parity are not degenerate in 
mass and there is an "almost" Goldstone boson with ~ rriq + 1/N. 



To be more specific, one can show |T2| that in QCD2{N oo) the 



behavior of the proper two-point correlation function is as follows: 

{0\T{qLqRix), gRgL(0)}|0) ~ x"^. (15) 

Such a behavior together with cluster property as x — * oo implies the ex- 
istence of the condensate at = oo in a full agreement with our previous 
discussion. At the same time, for any finite but large A^, the correlator falls 
off very slowly demonstrating the BKT-behavior with no signs of contradic- 
tion to the Coleman theorem. 

Having these general remarks on QCD2{N) in mind, we turn into our 
main subject. 

3 Large order behavior in QCD2{N = oo). 
3.1 't Hooft model. 

Let us consider the asymptotic limit = —q^ ^ oo of the two-point corre- 
lation function 0, [|T2| : 



I j t/xe^«^(0|T{g*75g(x),TO(0)}|0) = P{Q^). (16) 

It is clear, that the large behavior of P{Q^) is governed by the free, 
massless theory, where 

W-oo) = -^lnQ2. ^17) 
At the same time the dispersion relations state that 

pm = ^^^ E (18) 

and the sum is over states with even n because we are considering the pseu- 
doscalar currents. Here residues /„ are defined as follows 



{Q\q^^,q\n) = ^|^-^U, n = 0,2,4,... (19) 

Bearing in mind that for large n, and — > mQir'^n, we recover 

the asymptotic result ([171). We can reverse arguments by saying that in order 
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to reproduce InQ^ dependence in the dispersion relation ([TTD, the residues 
must approach the constant ~ {f^n+i — f^n) for large n. 
Now consider a Q~'^'' expansion for the correlator (|1^) in order to find 

the coefficients of this series at large k 

p(g^)~Ec2./^~Ec2.(^)^ (20) 

k k ^ 

As we mentioned earlier, in two dimensions the perturbative expansion Ck{g'^N)^ 
and the 1/{Q'^Y -expansion coincide. 

Now, if we knew /„ and m„ for arbitrary n we could calculate the sum (0) 
precisely, and thus, we would find the coefficients Ck from (|20|) . Unfortunately, 
we do not know them. However, the key observation is as follows: in spite 
of the fact that we do not know an analytical expression for /„ and m„ for 
arbitrary n we still can calculate the leading behavior of Ck- The reason for 
that is related to the fact that the only asymptotics of residues fn = ^, n — > 
oc and masses = mlTi'^n, — > cxo are essential; The corrections to /„ = 
1 + 0(l/?7,), = m^Ti'^n + 0(1), n ^ oo might change the preasymptotic 
behavior of c^, k oo, but can not change the factorial behavior (fc)!, found 
below. Using the asymptotic expressions for /„ and m„ we find that P{Q^) 
is expressed in terms of transcendental function "^(z) = where z = J^^^ ■ 
However we can trust only in the leading terms of the corresponding formula: 

\^ J ^ ^ 271 "^n + z n 27r 



lnz + 7^ + -^-E^-Tk 



2z 



k=^ 



2k 



^ _ (2fc)!2(- P^-^ 
moTT^ "^^ (27r 



232 ^2fc - (21) 



where B2k is asymptotic expression for the Bernoulli numbers. 

A few comments are in order. First of all, we have explicitly demonstrated 
that the coefficients C2k in the operator expansion J2c2k/Q'^'' are factorial 
divergent in high orders, C2k ~ (2A;)!, so the expansion is asymptotic in full 



agreement with the general arguments of ref. pT |. 

Additionally we note that the only even powers g'^/Q'^ are essential in the 
expansion ( generally speaking, arbitrary powers oi g^/Q^ could contribute). 
Nonleading terms in and might contribute to the odd powers g'^/Q'^. 

From the physical point of view this factorial behavior is related to highly 
excited states with excitation number 2nQ ~ Q'^/rriQ, and not with multiple 
production as one might naively expect. Indeed let us consider Q^/ml ^ z < 
2no in the expansion It is clear that the main contribution 

comes from k ~ 2no and (2?7,o)^"° ~ (2^o)! exactly corresponds with the 
behavior found above. It is in agreement with phenomenological analysis 
| TT[| , where it was assumed that the production of a highly excited resonance 



might be responsible for the large order behavior. 
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From the theoretical point of view we would expect that this behavior 
is related, somehow, to purely imaginary instantons (in order to provide the 
correct (—1)*^ behavior). An additional arguments in favor of this point will 
be discussed latter. 

We would also like to point out that the numerical coefficient, which 
enters to the formula (|2l|) is follows: 



C2kg ~ ^~Q2~) ■ ^^^) 

At the same time, the perturbative 2/c-loop graph gives a different contribu- 
tion ~ (^;^)^^ with an extra factor I/tt per each coupling constant. This 
extra factor tt must be taken very seriously as it is a large parameter. We 
definitely know (from exact solution), that the real scale of the problem is 
mpTT^, and not ttIq = g'^Nc/ir as one would naively expect from the pertur- 
bative theory. This means that vacuum condensates which are determined 
by non-perturbative physics come into the game. Even more, their contri- 
bution is much more important than pure perturbative diagrams. Let us 



note, that the lowest vacuum condensates, found exactly in [^, exhibit this 
additional factor vr. Thus, the factorial growth is related, somehow, to the 
nonperturbative physics. 

To further investigate the nonperturbative nature of the asymptotic series 
(in order to support the previous arguments) let us, instead of correlation 



function (|ly), consider the following difference of correlators: 



rfxe^'?^{(0|T{TO(x),TO(0)}|0) - (0|r{gg(x), gg(0)}|0)} = APiQ'). 

(23) 

One can argue that in the chiral limit rriq —>■ the perturbative contribution 
to (|23| ) iz zero. At the same time dispersion relations lead to the same result: 



the coefficients of the OPE are factorially divergent. This growth is related 
not to some perturbative diagrams, but to nonperturbative physics. We will 
present more arguments for this point of view in the next section. 

Finally we have explicitly demonstrated that the OPE is an asymptotic 
series. However we can not answer the important question of what kind of 
operators are responsible for such behavior. The reason for that is simple- 



too many operators contribute to the correlation function (|T^) and the cor- 
responding classification problem is quite involved. In the following we will 
consider a special heavy-light quark system, where such an identification can 
be made. We find that some vacuum condensates exhibit a factorial growth. 
Exactly this fact is the source of such an asymptotic behavior. 



8 



3.2 QCD coupled to adjoint fermions. Instantons. 

Now we repeat the preceding analysis for the much more interesting model 
of QCD with adjoint Majorana fermions pO| , ^, 

Sadj = J d^xTr[ - ^F^uF'"' + i^^^^D^^ + m^^] (24) 

As is known, the most important difference with t'Hooft model is that the 
bound states may contain, in general, any number of quanta. In other words, 
pair creation is not suppressed even in the large limit. The problem 
becomes more complicated, but much more interesting, because the pair 
creation imitates some physical gluon effects. 

We consider the following correlator analogous to ([T6|) : 

i f dxe*^^(0|T{-^TrM(x), -^Tr^'^ (0)}|0) = P2{Q^), (25) 

where ^ = \l/^7o and the label P2 shows the number of partons in the external 
source ^"^{x); the factor l/N^ is included in the definition of the external 
current in order to make the right hand side of the equation independent on 
A^. In the large limit the leading contribution to correlation functions is 
given as before by 

P2{Q^ ^ 00) = -^\nQ\ (26) 

The additional factor 2 comes from two options in calculation of Tr and 
related to Z2 symmetry mentioned in [Q. 

Now the problem arises. In t'Hooft model we definitely know that only 
2-particle bound states contribute to the corresponding correlation function. 
However, this is not true for the model under consideration and any states 
may contribute to P2. The key observation is as follows: any pair cre- 
ation (quantum loops which describe the virtual effects) is suppressed by 
a factor g'^Nc/ because of dimensionality of the coupling constant in two- 
dimensions (in a big contrast with real 4-dimensional QCD)^. Besides that, 
the quark mass term produces the analogous small factor m^/Q^ and can 
be neglected as well. Thus, information about highly excited states (which 
provides the logQ^ dependence) can be obtained exclusively from the anal- 
ysis of the correlation function at large Q^. In this case the analysis is very 
similar to 't Hooft case, considered in the previous section: 



1 - /(m^vr^) 
gineqnarray{0\—'^'^\ni) = \ rii ^ 1, n\tlZ 

Nr V VT 



■^Naively, one could interpret such a result that the mixing between different number of 
partons, is highly suppressed. Such a conclusion would be in contradiction with numerical 
results However, as we argued in recent papers |2|] the puzzle can be resolved by 

introducing a nonzero value for vacuum condensate (^'5). Such a condensation does not 
break any continuous symmetries. Thus, no Goldstone Boson appears as a consequence 
of the condensation. 
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= /„! = 1, rno = . (27) 

71 

The only difference is the doubhng of the strength of the interaction g"^ 2g'^, 
IP] and the additional degeneracy Z2, mentioned above. Now the formula for 
P2 (^Bp can be easily recovered: 

P2iQ^^oo)-pm = l E ^^r§--^MQ^), (28) 



ni=0,2,4... "^'1 + 



As before, any correction to the asymptotic expression (|^, such as /^^ = 
1 + 0(l/ni), m^^ = rngTr^ni + 0{mg) will produce power corrections ~ l/Q^ 
and they are not interesting at the moment. 



Apparently, the formulae (p6| - p8D are very similar to eqs. (|T7|-|l^) which 
correspond to the 't Hooft model. However, there is a big difference in 
interpretation of these two cases: In the 't Hooft model we have exclusively 
two- parton states (two bits, in terminology of refs.[pi,|20|]). They saturate 
the dispersion relations. 

In the case (p6| - p8D we have much more states with arbitrary number of 
partons. As we explained in [E^ the mixing between the different numbers 



of partons is not suppresed because of (\E'\E') -condensation. Effectively, how- 
ever, all these complex states contribute to the correlation function ( P^ ) in 
the same way as in the 't Hooft model. In this case the integer number ni 
from ( p7|) should be interpreted as an excitation number of 2- bits in those 
states. The matrix element /^^ ~ 1 can be interpreted as a total probabilty 
to find 2- bits among the complete set of the mixed states. The total number 
of states is increasing the mass increases. Thus, the probability to find 2 - 
bits in the given state is decreasing correspondingly. However, the dispersion 
relations (|28| ) tell us that the total probability /^^ with the given excitation 
number rii remains the same. 

We can repeat the previous analysis, which led us to the formula (pi|) 
with small changes. Instead of factor in the expression ( pT]) we will find 

B2k =^ 2^*''i?2fc 

for the theory with adjoint matter. This difference was mentioned above and 
is related to the doubling of the strength of the interaction g"^ —>■ 2g'^. 

How could one interpret this result? First of all, let us recall, that facto- 
rial behavior may occur for three different reasons: ultraviolet renormalons, 
infrared renormalons and instantons. Clearly, the first two reasons can not 
cause for such behavior in a field theory with the dimensional coupling con- 
stant g. Thus we expect that some kind of classical solution should be re- 
sponsible for such behavior. If we accept the instanton hypothesis, then from 
the very general arguments, one would expect that the instanton contribution 
with action S to the large fc-order coefficients Ck is given by |p3|| : 



Ck{g')' ikV.S'\gy (29) 
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In this case the factor 2, mentioned above has the following interpretation: 
when we go from the QCD with fundamental matter to the theory with 
adjoint matter the instanton action is decreased by a factor of 2. In this case 
the factor S~'' from the formula ( pUD is exactly equal 2^. It can be interpreted 
as the decreasing of instanton action by factor 2. Why the instanton with the 
action one-half is allowed in the theory with adjoint matter and forbidden 
in the theory with fundamental fermions ? This question has yet to be 
answered. 

Let us conclude this section by noting that from arguments given above 
we expect that some classical, pure imaginary solution (we call it instanton), 
is responsible for the factorial behavior found above. 



4 Heavy- light quark system in QCD2. 
4.1 General remarks. 

As we mentioned in previous sections we are not able to identify the n! be- 
havior in the OPE (found from the spectrum) with some specific vacuum con- 
densates. Such an identification can be done if one considers the heavy-light 
quark system. In this case the problem is reduced to the analysis of vacuum 
expectation value of the Wilson line (W) = {q{x)P exp{ig Jq Afj^dx^)q{0)). 
Indeed, if we consider as in [^,|]2S[ the correlation function {T{qQ{x), Qq{0)}), 



describing this system, we end up (in the limit Mq oo) with the object 
which is completely factorized (in accordance with HQET, see e.g. review 



26[| ) from the heavy quark: 

{T{qQ{x),Qq{0)}) ~ {q{x)P exp{ig I A^dx^)q{Qi)) + perturb, part. (30) 



By definition, (W) in this formula is understood as the Taylor expansion: 

ra=oo 1 

{W) = (0|g(x)Pe^^/o ^'"^--g(O)lO) = ^ 7^(^(^m^m)'"?) (31) 

n=0 l^^^ 

All nontrivial, large distance physics of the system is hidden there. Together 
with perturbative contributions one should expect the following behavior for 
this correlator : 



(T{gg(x),gg(0)})~e-^-^ (32) 

The perturbative terms, proportional to ~ (^^)"(a;^)", contribute to (^) as 
nonperturbative ones due to the dimensional coupling constant g in 2d, thus, 
they interfere with expansion (^TJ). As was suggested in the same context by 
Shifman one can get rid of the perturbative terms by considering a spe- 



cial combination of scalar and pseudoscalar correlation functions (analogous 
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to (|2^) with the replacement of a hght quark for a heavy quark). The per- 
turbative contribution vanishes in the chiral hmit for the such combination 
and we can study the pure nonperturbative physics. 

Let us note that the general connection (based on dispersion relations) 
between spectrum and vacuum condensates was considered earlier |]21]]. It was 
proven that the OPE is asymptotic series. Besides that, for demonstration 
purposes, it was suggested a specific model for the spectral density (linear 
trajectory) and it was found that the vacuum condensates (|3lD get the form: 
{q{x^D^f^q) ~ {2n)\. 

In this section we essentially follow the steps from the paper with the 



only difference being that we start from theory defined as QCD2 and derive 
the Qq spectrum from this Lagrangian. We find that we have not a linear 
spectrum but rather En ~ y/n in QCD2 for a heavy-light quark system. It 
gives a different behavior for the vacuum condensates (g(x^D^)^"g) ~ n\. 
However, the main statement that the OPE is an asymptotic series remains 
the same. 

Before going on, we would like to make the following remark. In QCD2 
one can calculate the appropriate vacuum condensates in the chiral limit from 



first principles |28]. Such a calculation (which will be reviewed in the next 
section) leads to puzzling results. Roughly speaking, the results of direct 
computation do not agree with indirect calculations based on the dispersion 
relations and spectrum. We formulate this puzzle as well as its resolution in 
the next section. Anticipating the event, we would like to note here that the 
origin of the puzzle is the factorial divergent coefficients in the asymptotic 
series. If these expansions were convergent series, we would expect an exact 
coincidence of the results, based on two these methods. However, the analysis 
in field theory demonstrates that this is not the case. 

4.2 Spectrum Qq system in QCD2{N = 00). 

As we have discussed in previous sections, if we knew the spectrum of highly 
excited states we would calculate (via dispersion relations ) the large order 
behavior for the corresponding correlation function. As we mentioned above, 
the heavy-light quark system is very special in this sense, because it allows us 
to identify the corresponding factorial behavior with specific vacuum conden- 
sates. This is the main motivation for the present section: find the spectrum 
for highly excited states. Let us note that the heavy-light quark system 
in this model was considered previously in ref.|29[, but in a quite different 
context. 

As we discussed in section 2, the spectrum of highly excited states in 
QCD2 is linear (^). This is certainly true, but only for the finite parameter 
mg, with n 00. We are now interested in a different limit, when rriQ —>■ 00 
first, and n ^ 1 afterwards. These limits do not commute and we have to 
start our analysis from exact original equation @. 
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In order to perform the limit ttlq ^ oo in 't Hooft equation let us 
make the following change of variables: 

X = 1 a. < a < — — = oo. (33) 

rriQ mo 

Additionally we would like to rescale the wave function and redefine the 
energy scale (from now on, the counting of the energy starts from mq) in the 
following way: 



ml = {niQ + c:^ tuq + ^niqEn, 0„(x) -> J— 0„(a). (34) 

y mo 

Once all these changes have been made we arrive to the following equation, 
which replaces the 't Hooft equation, for the heavy-light quark system in the 
limit mq —>■ oo. 

2E„0„(a) = mo a H ^ — 0n(a) - m^P / d/3- — , (35) 

a mo Jo [a — pY 

The new set of wave functions in terms of new variables is orthogonal and 
complete: 

V0n(a)0„(/3) = 5(a - /3), / rfa0„(a)0„,(a) = (36) 
n Jo 

For the future analysis we need not only the wave functions, but also some 
physical matrix elements in terms of these wave functions. It is convenient to 
separate the common factor, related to mq and define the matrix elements 
in the following way Q: 

/iV f°° — 

(0|gz75(5|n) = \ —^momqfn, fn= (ia0„(a). (37) 



Using the parity relation [Q, which in our notations takes the form, 

daMa) = 1}l r da^^ (38) 
mo "'0 a 

one can show that the scalar matrix elements have the same expression in 
terms of wave functions as the pseudoscalar ones (P?]). 

Having these results in mind, and using the standard technicques[|l|- 0, 
one can calculate |^the matrix elements fn and energies En in the quasi- 
classical approximation for n ^ 1: 

K = 2mov/^(l + 0(i^)), /„^ = ./^(l + 0(i^)). (39) 

n \ n n 



■*We keep the same notation /„ for the corresponding matrix elements. For hght quark 
system they are defined in the different way (0). We hope it will not confuse the reader. 
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This is the main result of this section. It will be used in what follows for the 
calculation of large order behavior and high dimensional condensates. 

Let us conclude this section with few remarks. First of all, as was ex- 
pected, in the uiq — > oo limit the equations PP] ), do not depend on mg 
(after an appropriate rescaling) in accordance with HQET (see e.g. review 

Our second remark is the observation that the chiral limit rriq = is very 
peculiar. In particular, one can not take the limit ruq = in the identity (|38D , 
because it clearly leads to a nonsensical result. The reason for that is very 
simple. The edge region a ~ rriq plays a very important role in making the 
theory (and this identity in particular) self consistent. As a consequence, one 
can not derive the boundary conditions on wave function from the truncated 
equation (|35| ) where the limit mq — >■ oo already has been taken. In order 
to do so we need to come back to the original 't Hooft eq.(0). We shall 
return to this point in the Conclusion. We believe that the situation in 
four dimensional QCD is quite similar in that the information about edge 
behavior in the theory can not be found from the truncated Lagrangian with 
the limit ttlq — > oo already taken. 

4.3 High order condensates from duality and disper- 
sion relations. 

The starting point, as usual, is the correlation function 

P(Q2) = , J e"''dx{0\T{qQ{x),Qq{0)}\0) (40) 



e''i''dx{0\q{x)Pe'^ -^0 ^^'^''^q{0)\0) + perturb, part. 



We follow and choose the external momentum go = ~ E, qi = 
very close to threshold. For positive values of E the correlation function can 
be written in the following way: 

poo A .1 

P(E)= e-^^dt{0\q(t)Pe Jo '^°'^^q(0)\0) + perturb, part. (41) 
Jo 

Let us note that we could consider the difference of two correlation functions 
(like (^3|)). In that case the perturbative contribution to ( Pl) would vanish. 
For large enough E ^ tuq one can expand eq. ( ^T]) in 1/E |^1|] : 

P{E) = ^[{qq) - -^{qPoq) + -^^{qP^) ' •••] ^ V^^rturb.part, (42) 

where Pq = ^-Dq is the time component of the momentum operator. 

Our goal now is to substitute the asymptotic expression (^) for the 
matrix elements and energies En into the dispersion relations (analogous 
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to (0)). These will determine the higher order corrections to the correlation 
function as well as the large n behavior of the vacuum condensates {qD'^'^q). 
The appropriate dispersion relation states that 

TV P N \ 
P(E) = —moV^y ^ -Y—-^ (43) 

where e is external energy measured in moi/vr units. In this formula we have 
taken into account the following key observation which we have used earlier in 
the deriviation of eg. (|2T|) : the corrections ~ 1/n to the asymptotic behavior 
of the residues /„ and energies En (^) might change the preasymptotic factor 
for the large order behavior. However, these corrections can not change the 
main result- the factorial growth of the coefficients found below. This is the 
reason why we can not calculate the corresponding coefficients exactly, but 
only the leading factorial factor. For the same reason we do not consider the 
special difference of two correlation functions (like (p3D), where perturbative 
contribution is exactly canceled. In this case if we knew /„, En exactly, we 
would calculate the nonperturbative condensates exactlyl Unfortunately, this 
is not the case. Thus, we ignore all complications related to the separation 
of pure nonperturbative contribution from the perturbative one. 

Let us note that the sum in eq.(^) is divergent at large n. This diver- 
gence is related to the necessity of a subtraction in the dispersion integral: 
P{E) —>■ P{E) — P(0). Besides that, at large energy E ^ mo one can es- 
timate the behavior P{E oo) ~ logE, which corresponds to the pure 
perturbative one loop diagram. These same features were present in our 
analogous previous formula (|2l|) . 

With these remarks in mind, our problem is reduced to the calculation of 
the coefficients Ck at large k in the following expansion: 

P(e) — subtractions ~ —= subtractions (44) 

-loge + ^Cfc^. 

k ^ 

We note that the only difference with the previous formula (^) is the de- 
pendence of the sum on ^/n rather than n itself. 

The way to evaluate the coefficients Ck is as follows. We are going to use 
the standard idea (see, e.g. [Q) to present the sum in terms of the integral: 

Y: F{n) = / -i^dx, Fix) = / f{x)c-^'dt (45) 
r^i -^0 e^-1 Jo 

However, in our case we have ^/n rather than n itself. The corresponding 
generalization is known as well: 

Y mn)] = r r ^^^^^P^dxdy, e-^^(P) = H h{x, y)c-^^dxm 
~^ Jo Jo — 1 Jo 
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For our particular case this formula gives the following representation for the 



sum 



n 



!/2 



Pie) ~ y , = --= I dx I dy^^^ (47) 



where we have introduced a constant c for the future convenience as an 
auxiliary parameter. Such a parameter is actually present in original formulae 
(^), however we believe that the main factorial dependence does not depend 
on it. Thus, after differentiating with respect to c we put c = at the very 
end of calculation. One more remark regarding formula (0). Only even n 
should be taken into account in this formula. However, by redefinition of 
parameters c and energy e, the problem can be reduced to the same integral. 
The result is an extra power dependence, which is beyound our scope of 
interest. Additionally, a subtraction which should be made in this formula 
to get a convergent result, has no influence on at large k (|4^) . 

The integration over x in the formula (0) can be executed using the 
following expansion: 

where are Bernoulli numbers. Bearing in mind that the appropriate 
integral from the (|47D is known exactly 



Jo 



-dx = \ -e~y^^, (49) 
X V c 



and replacing x'' from the formula (^) by (— 1)'^(^)'^, we arrive to the fol- 
lowing expression for the sum (^): 

PW^if ^ye-»..|§(-l)'(l)'(/ie-»v^). (50) 

In this formula we ignore few first terms (proportional to Bq,Bi) because 
they: a) do not contribute to large order coefficients c^. A; 1 , b) they are 
divergent and thus, require some subtractions discussed previously. 

The key observation is as follows: The nonperturbative part (which is 
our main interest) in the expansion (^2|) is determined by odd powers of 
l/(e)^"''^^. Such terms can be easily extracted from the formula (^) by 
expanding the exponent 

(51) 
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and executing the integration over y: 







?/' / -I- 1 

dye-^yy\ = (52) 



It is now clear, that only odd / = 2m — 1 terms in the formula ( |5TD contribute 
to the coefficients related to nonperturbative part (^21) . For small parameter 
c the coefficients can be easily calculated R by noting that for / = 2m — 1 

the appropriate terms from (|51| ) take the form ~ (2m-i)! • '^^^ nonzero 
result in this case ( after differentiating {-j^Y taking the limit c = ) 
comes exclusively from the term = m — 1 in eq.(^Ol). Finally we arrive at 
the following asymptotic expression for the odd coefficients Ck from the series 



1 

C2fc-i^- fc^2fc-i ' (53) 
Comparison with the original series (42) suggests that as dimension of the 



operator grows, their vacuum expectation values grow factorially, 

{qPo'l) - mi^^of'nl (54) 

From this formula one could naively think that half of the VEV's in (^) 
vanish because the odd Bernoulli numbers are zero. However, we think that 
this additional "selection rule" is accidental in its nature and thus, it should 
not be considered seriously. We believe that this accidental vanishing of 
the coefficients Ck in the OPE is related to our approximation ( ^9|) for the 
matrix elements /„ and energies En- Thus, we expect, that half of the VEV's 
which formally vanish in the leading approximation, are actually not zero, 
but suppressed by a factor ~ l/ra, where n is dimension of an opeartor. 

Having demonstrated the main result of this section, a few comments are 
in order. First of all, the OPE for P{e) is an asymptotic series as it must 
be in agreement with general arguments ||21]]. Besides that, we expect the 
same behavior for analogous vacuum condensates in four dimensional QCD^ 
27[] . Additionally, the scale for the vacuum condensates is mg. We shall see 



in the next section that different calculations of the same condensates give 
somewhat different results. We explain this puzzle later, but note for now 
that n\ behavior in (0) plays a crucial role in the explanation. 

Finally, as we shall see, the n\ behavior has very general, essentially kine- 
matical origin in the large limit. This property is related to so-called 
master field. 



^ As we mentioned before the factorial behavior does not depend on the particular 
magnitude c. However, for c — Q the calculations are much easier to present. Nonzero 
values of the parameter c might change the preassymptotic behavior, which is beyound 
this method. 
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High dimensional condensates from the the- 
ory. Puzzle. 



One can show [281, that the vacuum condensates which enter into the formula 



(^) at = oo in the chiral hmit niq — > in two dimensions can be reduced 
to a form which contains the field strength tensor ige^^G^u only. Indeed, 
the covariant derivatives (gD" ■ D^D^q) placed at the very right and at the 
very left (near the quark fields) can be transformed into the operator ige^uG^i, 
using the equation of motion D^'j^q = 0. To do the same thing with operators 
Dfj^ which is placed somewhere in the middle of the expression, we need to 
act, for example, on the right until the quark field is reached. By doing so, 
step by step, we create many additional terms which are either: commutator 
like [D^, D^] = —igG^i, which is the field strength operator or commutators 
like ~ [Dx,t^yG^y\. Fortunately, in two dimensions these terms are related 
to creation of the quark- antiquark fields and we discard them in the chiral 
limit because they do not give l/ruq enhancement, see formula ( [56| ) Q. We 
discuss the exact correspondence in the Appendix, but for now we emphasize 
the existence of factor n\ in the corresponding formula: 

{q{x,D^f^q) ^{x^YnmgEYq), n » 1, E ~ c.^Ga.. (55) 

Thus, we end up with the vacuum condensates {q{gEYq) which are ex- 
pressed exclusively in terms of the field strength tensor. Such vacuum con- 
densates can be calculated exactly in the chiral limit p8|. The reason for this 
incredible simplification is the observation that in QGD2 a gluon is not a 
physical degree of freedom, but rather is a constrained auxiliary field which 
can and should be expressed in terms of the quark fields. At the same time, 
in the large limit, the expectation value of a product of any invariant oper- 
ators reduces to their factorized values|3l|]. Exactly this feature of the large 



N limit (based on analysis of the so-called Master Field) makes it possible 
to calculate the vacuum condensates exactly. Our final expression takes the 
form p8| : 

^(g(^^7e,.G,.75)"g) = (-^)"(?"^) (56) 

We interpret this expression as follows: Each insertion of an additional factor 
proportional to field strength tensor gE, gives one and the same numerical 
factor (0). This situation can be interpreted as having a classical master 
field ||3l| which we insert in place of gE in the vacuum condensates. Because 
of its classical nature, it gives one and the same numerical factor. 

Secondary, it is important to note that the vacuum condensate of an 
arbitrary local operator can be reduced through the equation of motion and 
constraints to the fundamental quark condensate (|^). 

^ Of course this is not the case in four dimensions where [Z?^, G\a\ is an independent 
operator which can not be reduced to quark fields. 
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Finally, we must comment on the effective energy scale which enters into 
the expression (0): This is not ttiq as one could naively think, but rather 
mgjj = m^/iriq ^ rn^. The obvious technical reason for that is related 
to the fact in the light cone gauge (where the theory has been quantized) 
A_ = -^(^0 — ^i) = we have few constraints: the usual constraint in the 
gauge sector (Gauss law)[]: 

d^E'^'r^g{qi'ql-j^6'^''qir^), (57) 

with right moving fermions considered as dynamical degrees of freedom. 
The left-moving fermions g_ are non-dynamical degrees of freedom in this 
gauge. The latter can be eliminated by the following constraint: 

~ — g_ (58) 

(for more details see e.g. |1^). This relation explicitly explains the origin of 
the factor l/rriq which is present in the formula (^) (see for details). 

Before formulating the puzzle, we would like to pause here in order to 
explain the definition of the high dimensional vacuum condensates. As is 
known, they are perturbatively strongly ultra violet (UV) divergent objects. 

As usual, all vacuum condensates should be understood in a sense that 
perturbative part is subtracted. The subtraction is organized by introducing 
of the so-called normalization parameter fi. In general, vacuum condensates 
do depend on this parameter fi. The gluon condensates of dimensions four, 
six,eight,... in four dimensional QCD are perfect examples where perturba- 
tive parts are divergent, put nonperturbative parts ( the remnants, which are 
left after subtraction) are perfectly defined. One could naively think that the 
instantons might spoil this picture, because they give ultra -violet divergent 
contribution to high order condensates. However, we do not think this is 
the case and one can argue that the definition we have formulated remains 
untouched even with small size instanton effects taken into account. 

The argument is as follows: We should treat the small size instantons and 
the small size perturbative contributions on the same basis. So, we should 
subtract both these divergences at the same time in order to get the so-called 
" non-perturbative" condensates which define the large distance physics. Of 
course small-distances physics does not disappear when we do such a sub- 
traction. According to Wilson OPE the corresponding small distance contri- 
bution (perturbative part and small instanton contributions) should be taken 
into account separately. 

In two dimensions this problem of course is much simpler. However, the 
perturbative contributions to the condensates are divergent as well. Never- 
theless, the high dimensional condensates perfectly exist. Indeed, we have a 

''This formula explicitly demonstrates why quark condensate appears in ( ^6| ) each time 
when we insert an extra gluon field E"^''. 
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formula ([56|), with mixed condensates expressed in terms of chiral conden- 
sate {qq) where the latter are defined as the remnant after subtraction of the 
perturbative contribution (for the technical details, see |]16|). Thus, our for- 
mula (|5B|) is understood as a nonperturbative one, when we treat {qq) as the 



nonperturbative chiral condensate. Let us note, that the gluon condensate 
in this model is finite and can be calculated exactly [jl2| (see also the recent 



paper on this subject). In QCD2 it can be expressed in terms of the 
chiral condensate. 

With this remark in mind, we are now ready to formulate the puzzle. 
The scale which enters into the OPE (^2]) presumably is determined by the 
coefficients Ck from (|4^). The leading contribution to these coefficients can 
be expressed in terms of the vacuum condensates (^) 

~ (-^^)'(^?) ~ (ral/mgfkl, k » 1. 

3 

Thus, the characteristic scale of the problem is — . We see an extra factor 

' mq 

l/rriq in this scale. It has very clear origin ( pBD and might be very big in the 
chiral limit. At the same time, the characteristic scale which can be found 
from the spectrum ( PT| , ^2D , is much smaller and proportional to ttIq. This is 



the puzzle. 

The explanation of this apparent paradox is as follows. If our series 
(^,|^) were convergent ones, we would be in trouble, but fortunately, our 
series are asymptotic ones. Thus, in order to make sense to these series 
we have to define them and here we will use the standard Borell summing 
prescription f\ Once this prescription has been accepted, we can write 



down an expression which reproduces the asymptotic coefficients for large 
number k and at the same time is well defined everywhere: 



1 



We did not specify the parameter 7 in this equation on purpose. Suppose, 
we have a large scale in the problem determined by the vacuum condensates 
(P^). In this case dimensionless parameter 7 has a large factor l/m^: 

7 (60) 

The prescription (^) in this case states that the sum of leading terms ~ 
{l/iriq)^ gives a zero contribution (in the chiral limit) 

Ecfe7'-/ e-7-.0(atm,^0), 7-— (61) 

V J 7 (7 + 7 ) 



^We already mentioned in Introduction that the Borcl summabihty or its loss is not the 
crucial issue However, for the sake of definiteness (and for simplification) we assume 
in general discussion which follows, that series is Borel summable 
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to the physical correlation function! We would like to note that the effect 
does not crucially depend on the factorization properties for the conden- 
sates (|56|) as neither on our assumption of exact factorial dependence of the 
coefficients = k\{—)'' (^). Both of these effects presumably lead (apart 
to k\) to some mild fc-dependence which can be easily implemented into the 
formula (^) by introducing some smooth function /(7) whose moments 

Cfe = ^/ /(7)7-'-'exp(-i)d7 - 1 

exactly reproduce a k- dependence of the coefficients as well as of the con- 
densates. If this function is mild enough, it will not destroy the relation (|6T1) , 
but might change some numerical coefficients. 

If we could calculate the vacuum condensates (^) exactly (and not only 
the leading terms at rriq — > 0), we would find that the terms of order one give 
contribution of order one to the correlation function. Thus, subleading terms 
play much more important role in the final formula than the leading ones ~ 
l/rriq. The origin of this mystery of course is the factorial growth coefficients 
in the series. This observation actually resolves the puzzle announced in the 
beginning of the section. 

We would like to make a few more remarks regarding this subject in 
Conclusion which follows. 

6 Conclusion. 

We conclude this manuscript with the following lessons. 

a) . The OPE is asymptotic series. 

b) .We were lucky in a sense that in QCD2 the scales for high dimensional 
condensates ~ and for the spectrum ~ 1 were parametrically different. 
This was the reason why we noticed the difference very easily. The lesson 
we can learn from this example is that numerically leading terms in the 
asymptotic expansion may give somewhat negligible contribution into the 
final expression. We believe that this is not a specific result for 2d physics, but 
rather is a general property of a field theory associated with an asymptotic 
series. 

c) .Of course, we do not expect that in the real QCD^^ the vacuum conden- 
sates might exhibit some parameters similar to ~ 'ti"^. However, it might 
happen that some subseries of condensates possess a large numerical factor, 
let us call it L ^ 1. In this case we would expect that the corresponding 
contribution into the final formula will be suppressed, L~^. At the same time, 
summing of a subseries of terms which are order of one, would lead to the 
result which is order of one. This observation raises the following question. 

d) . Suppose we have a condensate (O^) of dimension k which has both 
parts: the enhanced part, proportional to and the "regular" one of order 
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1*^. As we have learned, the subseries which has big factor hke do not 
contribute much to the final expression. At the same time they are the main 
contributions into the condensate with the given dimension k. 

The moral is: when we use truncated expansions or approximate ap- 
proaches (like QCD sum rules), we inevitably study not the actual conden- 
sates, but rather, some effective condensates |^. One simple consequence 
of this is that the lattice calculation of vacuum condensate might be different 
from QCD sum rules analysis. 

e) .One may wonder what is the role of the scale l/rrig in this model? 
The answer is very peculiar. Indeed this scale can not be seen in spectrum 
however, from exact identities like ( |38D one can see that the edge region 
of order plays a very important role in maintaining selfconsistency of 
the theory. The calculation of the condensate where the region x ~ iriq 
gives whole answer is another example of the same kind. Exactly this infrared 
region determines the scale for the condensates (|56D , but not the scale for 
integral characteristics like the spectrum itself. See also remark after (^) on 
the same subject. 

f) . Since OPE is asymptotic series, it is good idea to keep only a few first 
terms in the expansion (like people do in the standard QCD sum rules ap- 
proach, |]53[) and to stop at some pointQ.Any hopes to improve the standard 



QCD sum rules (like the idea advocated in |Q) by summing up of the certain 
subset of the power corrections, and ignoring all the rest, might lead to the 
results which are much worse than the ones which follow from pre-improved 
version of the approach. At least in QCD2 such a procedure, as we learned, 
gives a parametrically incorrect result. 



We would like to make two more remarks which were not our main sub- 
jects, but look like interesting byproducts worthy of mention. 

g) .We argued that the n\ behavior found in QCD2 is related to instan- 
tons. Even more, the action of the instanton in the theory with adjoint 
matter is a factor of two less than in theory with fundamental quarks. The 
explicit realization of such a solution is still lacking. 

h) .We analyzed the Qq system in QCD2. We found that the edge behavior 
of the system is very peculiar and can not be found from the truncated theory, 
where limit mq = 00 is already taken. We believe that the same behavior is 
an inevitable feature of the real QCD4. 

thank Michael Peskin for the discussion of this subject. Actually, he raised the 
question of the effective nature for the condensates before this work was presented at the 
SLAC theory seminar. 

^°In asymptotic series with coupling constant A one should stop with the number of 
terms of order A~^. In particular, in case of QCD sum rules, where A ~ 1/3-;- 1/5 to 
be determined by the scale where power corections are 20 30%, one can estimate the 
maximum number of terms in the expansion is about 3 5. 
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8 Appendix. 

The main goal of this Appendix is a demonstration of the factor n\ in formula 
(|6^) in the chiral limit rriq = 0. 

{q{x,D,f''q) r^{x^rn\{q{gErq), n » 1, E ex.Gx.. (62) 

We sketch the idea only, ignoring for simplicity all normalization factors. 
We use the standard representation for 7^ matrices satisfying 

{lt^^,y} = 25-/.^, 7o = 7i = -i(^y, 75 = c^s, (63) 
7± = 7o±7i, D± = Do±Di, 7^ = 0, 7+ = 0, 7±7t ~ 1 =F 75 

Dirac equations take the following form in the chiral limit 

(7+/^-+7-^+)g = 0. (64) 

By multiplying these equations by 7-1- , we arrive to the following relations 

(l±75)D±g = 0, gD^(l±75)= 0. (65) 

Let us present our original condensate in the following form 

{qix,D,rq) ~ (g(i±^ + i^)(x+D_ + x^D^y^q). (66) 

From Lorentz invariance it is clear that a nonzero result ~ (x_x+)" comes 
from terms with equal numbers of D_ and with all possible permutations. 
Our problem is how to count them. First, we pick up the projector (1 + 75) 
from the expression (|66D . According to (|65|) we should move to the right 
and to the left (to reach the quark field) using the following commutation 
relations 

[D+,D_]~e^,G'^, (67) 

Here the field strength E can be considered as constant (not operator), be- 
cause acting the operators D± on E leads to the pair creation. We neglect 
everywhere such terms in the limit — > 00, mg —>■ 0. The contribution with 
projection (1 — 75) gives exactly the same result after relabelling D_ 
and repeating the procedure described above. 
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Now one can see that these calculations are very similar (in the algebraic 
sense) to the oscillator problem with ladder operators satisfying the standard 
relations 

[a,a+] = l, a|0) = 0, (0|a+ = 0. (68) 

Our problem is the calculation of mean value of the operator x^" ~ {a + a^Y^ 
for the ground state: 

(0|x2"|0) ~ (0|(a + a+)2'^|0) ~ r(n + ^) ~ n!, n > 1. (69) 



This concludes our explanation of the factor n\ (|55D which we heavily used 
in our previous discussions. 



Let us repeat again that all factors from (|55| , |56D have clear meaning 
and can be explained without detailed calculations. Three steps are in- 
volved: l.The transition from the operator (g(x^D^)^"g) to the operator 
{q{gE)"'q) with a factor n\ was explained in this appendix. 2. The idea of 
Master Field predicts that each insertion of E into the expression for the 
operator {q{gE)"-q) gives one and the same constant M^.r. 3. The constraints 



(0,|^) explicitly demonstrate that M^j^j ~ ^ 
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